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Abstract 

We evaluate the fluxes of up-going muons detectable in a neutrino telescope and due 
to the annihilation of relic neutralinos in the Earth and in the Sun, taking realistically 
into account the fact that neutralinos might provide only a fraction of the local (solar 
neighbourhood) dark matter (DM). We determine the relation between the exposure At 
of the neutrino telescope {A being the telescope area, t the live time) and the explorable 
range of the neutralino mass and composition, when the signal-to-background ratio is 
optimized by an appropriate angular selection. 
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1 Introduction 

The perspectives of observing, with neutrino telescopes of large area, indirect signals 
of relic neutralinos form celestial bodies (Earth and Sun) are very appealing and are 
currently the subject of much investigation [1-3]. The sequence of the physical processes 
which could produce these signals are: i) capture of relic neutralinos by the macroscopic 
bodies; ii) subsequent accumulation of these particles in a region around the centre of these 
celestial objects; iii) pair annihilation of the accumulated neutralinos which would generate, 
by decay of the particles produced in the various annihilation final states, an output of 
high-energy neutrinos. These neutrinos may be detected in a neutrino telescope in two 
ways: either by looking at events due to the 2v-interactions inside the detector (contained 
events) or by measuring the muons which are produced by z^^ and z/^ interactions in the 
rock around the detector and then traverse it upwardly. In the present paper we only 
discuss the latter case; we address the problem of the theoretical predictions for the fluxes 
of these up-going muons in a neutrino telescope and of their discrimination against the 
background due to atmospheric neutrinos. The present analysis aims at improving the 
previous ones in many instances: a) The calculations of the predicted signals for both 
macroscopic bodies (Earth and Sun) are presented in a common theoretical framework; 
in particular, this enables a comparative discussion about the informations obtainable by 
measurements of the two fluxes, b) It is realistically taken into account the fact that 
neutralinos might provide only a fraction of the local (solar neighbourhood) dark matter 
(DM) density. In order to establish whether and by how much the total DM density has to 
be rescaled, the neutralino relic abundance is explicitly evaluated at any point of the model 
parameter space where signals are calculated, c) It is flnally discussed the relation between 
the exposure At of a neutrino telescope (A being the telescope area, t the live time) and 
the explorable range of the neutralino mass (and composition), taking into account the 
optimization of the signal-to-background ratio (S/B). 

The theoretical framework adopted in the present paper is the one provided by the 
Minimal Supersymmetric extension of the Standard Model (MSSM) [4], which, imple- 
mented by a general grand uniflcation relation, may be formulated (for the neutralino 
sector) in terms of only three parameters: two mass parameters M2 and /i, and a third 
parameter tan/? defined as tan/3 = Vu/vd {vu and Vd are the v.e.v.'s that give masses 
to the up-type and down-type quarks, respectively). The neutralino, is defined to be 
the lowest-mass linear combination of the two gauginos (photino and zino) and the two 
higgsinos 

X = ai7 + a2Z + a^H^ + 04^2 • (1) 

Here 7 and Z are the fields obtained from the original flelds for the U(l) and SU(2) neutral 
gauginos, B and W3, by a rotation in terms of the Weinberg angle. 

It is useful to introduce a composition parameter P = ai^ + which gives the 
gaugino fractional weight. The neutralino mass m^^ and the coefficients a^'s depend on the 
three model parameters M2, jj, and tan/?. Since almost everywhere in the parameter space 
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there is a one-to-one correspondence: {M2, ji) <-> (Pjin^), at fixed tan/3, in the following 
we will use as a set of independent parameters: P, and tan f3. This allows us to present 
our results in terms of two quantities, P and m^, which have a direct physical meaning 
for the neutralino. In the following, for tan/?, unless differently specified, we will use the 
representative value tan/3 = 8. 

We still have to state which values we assign to the unknown masses of the Higgs 
bosons and of the sfermions; in fact these masses enter in the theoretical calculations, since 
Higgses and sfermions play a crucial role in the physical processes initiated by neutralinos. 
As for the neutral Higgs bosons, we recall that in the MSSM there are three neutral Higgs 
particles: two CP-even bosons h and H (of masses nih, mu with mu > mh) and a CP- 
odd one A (of mass niA)- Once a value for one of these masses (say, mh) is assigned, 
the other two masses (m^, m-H) are derived through mass relationships depending on 
radiative effects. In the present paper the Higgs mass m/j is set at the value of 50 GeV. 
As for the sfermion masses, these are set here at the value "7.^ = 1.2 m^, when > 45 
GeV, m J = 45 GeV otherwise. Only the mass of the stop quark is assigned a larger value 
of 1 TcV. The top mass is fixed at mt — 170 GeV. 

2 Neutralino local density 

An important quantity that enters in the capture rate of relic neutralinos by macro- 
scopic bodies is the neutralino local density p^. One can assume that is equal to the 
local value of the total DM density pi, when the neutralino relic abundance (O^/i^) turns 
out to be at the level of an (n/i^)min consistent with pi. On the contrary when {Q^h'^) 
is smaller than (n/i^)inin the value to be assigned to p^ has to be appropriately reduced. 
Thus we evaluate Q-^^h^ and we determine by adopting a standard rescaling procedure 
[5]: 
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Px~Pii when Q,^h > {flh ) 

Px = Pi (^u2\ — ' '^^^^^ VL^h? < {VLh'^)^iry 



Here (ri/i^)min is set equal to 0.03. For p;, we have used the value pi = 0.3 GeV cm ^. 

The neutralino relic abundance Vl^h? is evaluated here following the procedure illus- 
trated in Ref. [6]. By way of example in Fig.l we display fi^/i^ as a function of m^^ for three 
representative neutralino compositions: i) a gaugino-dominated composition [P — 0.9), ii) 
a composition of maximal gaugino-higgsino mixing (P = 0.5), iii) a higgsino-dominated 
composition {P = 0.1). As expected, out of the three compositions displayed in Fig.l, the 
gaugino-dominated state provides the largest values of Vt^h"^ . More substantial values of 
Vl^h? occurs if one considers purer gaugino compositions (P > 0.99), as shown in Fig. 2, 
where, together with the reference value P = 0.5, are also displayed the rather extreme 
cases: P = 0.01 (very pure higgsino composition) and P = 0.99 (very pure gaugino com- 
position). The very pronounced dips in the plots of Fig.s 1-2 are due to the s-poles in the 
X-X annihilation cross section (exchange of the Z and of the neutral Higgs bosons) . 
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Displayed in Fig.3 are the values of fi^/i^ versus in the form of scatter plots 
obtained by varying M2 and over a grid of constant spacing in a log-log plane over the 
ranges 20 GeV < M 2 < 6 TeV, 20 GeV < < 3 TeV. By comparing Fig.3 with Fig. 2 one 
sees that the large spread in values for Q^h^, displayed in Fig.3, is due to configurations 
of extremely pure composition. 



3 Capture rates and annihilation rates 

The capture rate C of the relic neutralinos by a macroscopic body may be evaluated 
by the standard formula [7] 



c = ^J2^^(MBmvi^hx„ (3) 

where is the neutralino mean velocity, a"ei,i is the cross section of the neutralino elastic 
scattering off the nucleus i of mass m^, Msfi is the total mass of the element i in the 
body of mass Mb, ('t'esc)i square escape velocity averaged over the distribution of 

the element i, is a factor which takes account of kinematical properties occurring in 
the neutralino-nucleus interactions. The elastic x-nucleus cross sections are evaluated 
according to the method presented in Ref. [8]. 

The annihilation rate Ta is expressed in terms of the capture rate by the formula [9] 



r^ = ^tanh2(^) (4) 

where t is the age of the macroscopic body {t = 4.5 Gyr for Sun and Earth), ta = 
(CC^)~^/^, and Ca is the annihilation rate per effective volume of the body, given by 



^ < av > , niy ,0/0 

Ca = ( — )^ (5) 

^ Vo 40GeV^ 

Vo is defined as Vq = i'^i^^'piT/ (2px 10 GeV))"^/^ where T and p are the central temperature 
and the central density of the celestial body. For the Earth (T = 6000 K, p = 13 g • cm~'^) 
Vo = 2.3xl025cm3, for the Sun (T = 1.4xlO^K, p = 150g-cm-3) Vq = e.GxlO^Scm^. cr is 
the neutralino-neutralino annihilation cross section and v is the relative velocity. < av > 
is calculated with all the contributions at the tree level as in Ref. [6], with the further 
inclusion here of the two gluon annihilation final state [10]. 

We recall that, according to Eq.(4), in a given macroscopic body the equilibrium 
between capture and annihilation {i.e. Fa ~ C/2 ) is established only when t > ta- It is 
worth noticing that the neutralino density p^, evaluated according to Eq.(2), enters not 
only in C but also in ta (through C). Therefore the use of a correct value for p^ (rescaled 
according to Eq.(2), when necessary) is important also in determining whether or not the 
equilibrium is already set in a macroscopic body. 
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In Fig. 4 we give the results of our calculations for tanli^(t/r>i) and for Ta in the case 
of the Earth for three representative x-compositions: P = 0.1,0.5,0.9. For simplicity, in 
these figures, as well as in the following ones, only the results concerning positive values 
of the parameter /j, are shown; similar results hold for negative values of The upper 
part of Fig. 4 shows that the equilibrium between capture and annihilation is not reached 
for m-)^ > mw'i thus a substantial suppression occurs in at large values of m-^ because 
of the factor tanh^(t/T^). In Fig. 5 the two quantities tanh^(t/Tyi) and F^ are plotted for 
other representative values of neutralino compositions: P = 0.01, 0.99 (the value P = 0.5 
is shown again for comparison). In Fig.s 4-5 we see that, for < 70 GeV, F^ shows 
the characteristic bumps due to fact that the capture is very efficient when the neutralino 
mass matches the nuclear mass of some important chemical component of the Earth (O, 
Si, Mg, Fe). Superimposed to these properties are the features due to the rescaling of p^. 
In fact this rescaling implies here both a general suppression and the appearance of the 
peculiar dips reminiscent of the singularities in the annihilation cross section (discussed in 
Sect. 2). For comparison with Fig. 4 we show in Fig.6 how the two quantities tanh^(t/ryi) 
and F^ would appear, if the rescaling were not applied. 

Whereas for the Earth the equilibrium condition depends sensitively on the values of 
the model parameters, in the case of the Sun equilibrium between capture and annihilation 
is reached for the whole range of m^, due to the much more efficient capture rate implied 
by the stronger gravitational field [7]. Fa for the Sun is shown in Fig. 7. 

What we have defined above is the annihilation rate referred to a macroscopic body 
as a whole. This is certainly enough for the Sun which appears to us as a point source. On 
the contrary, in the case of the Earth, one has also to define an annihilation rate referred 
to a unit volume at point r from the Earth center 

^A{r) = ^<av> n\r) (6) 
where n{r) is the neutralino spatial density which may be written as [9] 

n(r) = noe-'*^^'''; (7) 
here a — 2ttGp/ (ST) and no is such that 

Ta=^ <(7V > I cZVn^(r) . (8) 



4 Neutrino fluxes 

The differential neutrino flux due to x~X annihilation in a distant source (Sun) is 
given by 
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Z^^x/^ (9) 



dE^ And^ ^ dE^ 

F,f 

( F) 

where d is the distance from the source, F denotes the annihilation final states, 
denotes the branching ratios into heavy quarks, r lepton and gluons in the channel F; 
dNfjj/dEjy is the difi'erential distribution of the neutrinos generated by the hadronization 
of quarks and gluons and the subsequent hadronic semileptonic decays. 

In the case of the Earth we define a differential neutrino fiux generated at a distance 
R by the x^X annihilations which take place at a point r from the center of the Earth 

dN^ _ <av> ■n?{r) s.- ^{f) dNf^ 
dE^d^r ~ 87ri?2 ^ ^E^ (10) 

This is the neutrino spectrum which is used in our Monte Carlo code to evaluate the flux 
of up-going muons, as illustrated in Sect. 5. 

However, it is also useful to write the neutrino differential flux at a point on the Earth 
surface, once the integration over the distribution n{r) is performed, at fixed nadir angle 
^[7] 

d'^N. _ Ta sr- iF)dNf. 

dE^dcose ~ ^ UttRI ^^^^ dE, ■ (11) 

G{9) is given by 

G{9) ~ 2{2mp)e-'^""^ ^ (12) 

where (3 = ai?^, for > 10 GeV and 6 < 60°. This flux, as all other fliixes discussed 
afterwards, refer to surfaces orthogonal to the incoming direction. 

As far as the differential neutrino spectra are concerned, we briefly quote the relevant 
ingredients that we used in the calculation. A Jetset 7.2 Monte Carlo [11] code was used 
to compute the neutrino spectra due to b and c quarks, tau lepton and gluons. Following 
Ref. [12] we have neglected the contributions of the light quarks directly produced in the 
annihilation process or in the hadronization of heavy quarks and gluons, because these 
light particles stop inside the medium (Sun or Earth) before their decay. For the case of 
the Sun also the energy loss of the heavy hadrons in the solar medium was considered. 
The spectra due to heavier flnal states, i.e. Higgs bosons, gauge bosons and t quark, were 
computed analytically by following the decay chain down to the production of a b quark, c 
quark or a tau lepton; the result of the Monte Carlo was used to obtain the flnal neutrino 
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output. Because of the high column density of the solar medium, the absorbtion and the 
energy loss of the produced neutrinos were also included. 

As an example, we show in Fig.8 the neutrino spectra for neutralino annihilation in the 
Earth, for P — 0.5 and for three different values of the neutralino mass: = 20, 40 and 
120 GeV. The different shapes depend on which final states dominate in the annihilation 
cross section. For instance the spectrum at = 120 GeV shows the typical features due 
to the opening of the final state in the x^X annihilation. The relative magnitudes 

of the spectra at different values of are driven by the dependence of the annihilation 
rate on the model parameters. 

5 Signals of up— going muons 

Let us turn now to the calculation of the various distributions and rates for the up- 
going muons generated by the neutrino fluxes of Eq.s(9-10) (signals) and by the neutrinos 
produced in the Earth atmosphere (background). The conversion of the muon neutrinos 
into muons and the subsequent muon propagation up to the detector have been treated 
using a Monte Carlo simulation. 

The muon flux at the detector depends on the following quantities: the neutrino energy 
spectrum and spatial distribution, the differential v^ — N charged current cross section for 
the muon production and the energy losses and the multiple scattering of the muon in the 
rock surrounding the detector. 

As for the signals, the neutrino energies have been generated in the range 2 GeV 
< El, < according to the distributions of Eq.(9) for the Sun and of Eq.(lO) for the 
Earth. We emphasize that in the latter case the neutrino production point has been 
generated in a region aroimd the Earth center with a gaussian distribution ?i^(r), where 
n(r) is given in Eq.(7). The atmospheric neutrinos have been generated upwardly according 
to the differential distribution given by [13] for neutrinos with E,^ < 100 GeV and by [14] 
for neutrinos with E^, > 100 GeV. 

In the cross section for the charged current — N deep inelastic scattering we have 
employed the CTEQ parton distribution functions [15]. The energy losses and multiple 
scattering of the muon have been calculated by propagating the muon in standard rock by 
steps of 10^ g cm~^ when the muon energy was larger than 30 GeV, of 500 g cm~^ when 
it was between 1.2 and 30 GeV, and equal to E^/ (2.4-10""^ GeV) g cm~^ for the flnal step 
when < 1.2 GeV. The residual energy and the muon direction were kept in memory 
for each step from the muon production to the muon stopping. The detector entry point 
has been randomly generated with a uniform probability along the whole muon track. The 
muon energy and direction at the entry point has been calculated by interpolating between 
the initial and flnal values for the step where the entry point occurred. A random angle 
has been included to take into account a detector angular resolution of 2°. 



6 



6 Results and conclusions 

We show in Fig.s 9-11 a sample of our results for the signals expected from the Earth at 
the usual representative point tan = 8 (for the fliixes and the angular distributions shown 
in these figures and in the following ones the threshold for the muon energy is set at the 
value E^^ = 2 GeV). Fig. 9 displays the integral muon flux versus for various neutralino 
compositions. This has been obtained with a no-straggling approximation [16], i.e. the 
muons are produced and propagated in the forward direction (when this approximation 
is applied, then Eq.(ll) for the neutrino distribution is useful). The general structure of 
the flux clearly reflects the features of the annihilation rate in Earth, previously discussed. 
The horizontal line represents the Kamiokande upper bound for signals coming from the 
Earth [2]: 4.0 • IQ-^^cm-^s-i (90 % C.L.); this value refers to a half-aperture of Ot = 30° 
for a cone centered at the nadir. By comparing this upper limit with our fluxes we see that 
the regions explored by Kamiokande (at our representative point: tan P = 8,m/i and 
chosen as stated in Sect. 1) concern the neutralino mass range 50 GeV < < 80 GeV, 
for 0.1 < P < 0.9; the mass range is narrower for purer neutralino compositions. 

By way of example we show in Fig.s 10-11 the differential spectra and the angular 
distributions for the case P = 0.5 and = 120 GeV. In Fig. 10 the spectra obtained 
with the MC code for the expected signal and for the background are displayed by the 
two histograms. The solid line represents the spectrum for the signal obtained with the 
no-straggling approximation; it is clear from the plot that this approximation is quite good 
for the spectrum, and then also for the integrated flux. The spectra shown in this figure 
refer to the muons entering the detector within a cone of half-aperture 9t = 30° centered 
at the nadir. In Fig. 11 the angular distribution due to the signal, markedly peaked around 
the nadir direction, is compared with the almost flat background distribution. 

Now a very relevant question is what is the minimal exposure At necessary to provide 
a good S/B discrimination. To put it into a more quantitative way, we define an {At)min as 
the minimal exposure necessary to see a 4cr effect (with a signal of at least 4 events) [17]. 
(^t)min is displayed in Fig. 12a (the angle 9t of half-aperture from the nadir is again set 
at the value 9t — 30° here). However, it is important to remark that, due to the different 
features of the angular distributions for the signal and for the background (illustrated 
for instance in Fig. 11), (At) min may be optimized by appropriately choosing the angular 
window, i.e. Of. We call {At)opt the least value of (^t)min when 0^ is varied, and we denote 
by 6'opt the corresponding value of 9t. This optimization procedure has to be applied to 
{At)niin ior each value of the integral flux For the case illustrated in Fig. 11 we find 
^opt — 5°, (At) opt ~ 3020 m^ yr. Fig. 12b displays the dependence of (^t)opt on and P 
for our usual representative case. To obtain the result of Fig. 12b the optimization, based 
on the angular distribution, has been performed by using the results of the MG simulations 
at 4 values of the neutralino mass: m^^ = 20, 40, 60, 120 GeV, by interpolating between 
these values and by conservatively using for > 120 GcV the same angular window 
as for m^^ = 120 GeV. For the integral flux we employed the results (shown in Fig. 9a) 
obtained with the no-straggling approximation. By comparing Fig.s 12a and 12b one 
sees that the optimization improves signiflcantly the experimental capabilities especially 
at large values; in fact for these values the extension of the annihilation region is very 
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much concentrated around the Earth center. Fig. 12b shows reahstically by which amount 

the explorable ranges of the two neutrahno parameters and P increase as the exposure 
At is increased. We remind that the largest value for At obtained up to now is the one of 
Kamioka: At — 770 m^ yr [2] . A sizeable improvement in At will be provided by MACRO 
[18]. 

Let us turn now to the signals from the Sun. In Fig. 13 the relevant flux of up-going 
muons is plotted versus m^^ in the case of the representative point under discussion; $p 
is the flux through a surface orthogonal to the Sun direction. The horizontal line denotes 
the present upper bound of 6.6 • 10~^''cm~^s~^ (90 % C.L.) (an half-aperture 9 = 25° was 
taken here for a cone centered at the Sun direction) [2]. It is clear that this limit does 
not introduce any constraint here. In Fig. 14 the corresponding {eAt)-m-)^ plot is displayed 
(a cone of 5° half-aperture around the Sun direction has been used for the evaluation of 
the background, e is the on-source duty factor). In this plot we see that, in order to start 
some exploration about relic neutralinos using the flirx from the Sun, one needs at least 
(eAt) ~ 10^ m^ yr (this could provide information in the range 70 GeV< < 150 GeV). 
However, it is remarkable that by slightly increasing (eAt) above this level, the explorable 
range in terms of P and expands quite rapidly towards the high values of the neutralino 
mass. 

An interesting point concerns the relative importance of the up-going muon fluxes 
from the Earth and from the Sun. As we see from Fig.s (9,13), for light neutralinos the 
signal from the Sun is much smaller than the signal from the Earth; as increases the 
signal from the Sun may overcome the other, since the Sun gravitational field is much more 
efficient in capturing the neutralinos. More relevant for a comparative discussion about 
the signals from the two sources are the two diagrams of the exposure versus m-^ displayed 
in Fig.s 12b, 14. A very appealing situation is the one where the signals from both sources 
could be detected. For instance, assuming e ~ 0.5, an exposure At ~ 10^ m^ yr would be 
sensitive to both signals in the range 35 GeV < m,^ < 200 GeV. 

A word of warning is in order here to remind that the present evaluations depend 
sensitively on a number of free parameters. As far as tan/3 is concerned, we may notice 
that the dependence of the signals on this parameter is rather involved, but, as a general 
trend, the expected signals increase with the value of tan/3. For instance, for the Sun, 
moving from the value tan/3 = 8 to the value tan/3 = 20 the flux is enhanced by about 
a factor 5 for > 200 GeV; but, if the value of tan/3 is reduced from tan/3 = 8 to 
tan /3 = 2, the flux is suppressed (for > 50 GeV) roughly by an order of magnitude. The 
quantities calculated in the present paper are also rather sensitive on the values employed 
for the Higgs and the sfermion masses. For instance, our signals for the Earth are roughly 
proportional to m^^. In the present paper we have used for the Higgs and the sfermion 
masses values close to their present lower bounds. This choice enables one to establish, 
by use of the figiires previously reported, what is the minimum exposure necessary for 
a detector to be sensitive to the Earth and/or Sun signals and how the exploration into 
the DM neutralino properties may be progressively improved by expanding the size of the 
experimental apparatus. 
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To summarize we may conclude that the exploration in the realm of the relic neu- 
tralinos by neutrino telescopes, after the first insights provided by Kamiokande [2] about 
the Earth signals, will certainly be significantly developed by the underwater (under-ice) 
detectors under construction [3]. In order to have good perspectives for a systematic in- 
vestigation of this field one should finally aim at an apparatus of very large area of order 
of 10^ m^. 
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Figure Captions 

Figure 1 Neutralino relic abundance fl^h"^ as a function of for three different 
neutralino compositions P = 0.1 (dotted line), 0.5 (short-dashed line), 0.9 (long-dashed 
line). The upper plot refers to positive values of the lower to negative values. 

Figure 2 Same as in Fig.l, except that now composition are: P = 0.01 (dotted line), 
0.5 (short-dashed line), 0.99 (long-dashed line). 

Figure 3 Scatter plot for the neutralino relic abundance fl-^h^ as a function of m^. 
The two MSSM mass parameters M2, fJ, are varied in the ranges: 20 GeV < \iJ,\ < 3 TeV, 
20 GeV < M2 < 6 TeV. 

Figure 4 tanh^(t/TA) and for the Earth are given as functions of m^^, for the three 
representative neutralino compositions P = 0.1 (dotted line), 0.5 (short-dashed line), 0.9 
(long-dashed line). 

Figure 5 Same as in Fig. 4, except that now composition are: P = 0.01 (dotted line), 
0.5 (short-dashed line), 0.99 (long-dashed line). 

Figure 6 Same as in Fig. 4, except that now the neutralino local density p-^ is not 
rescaled. 

Figure 7 Annihilation rate F^ for the Sun as a functions of m-^, for the three rep- 
resentative neutralino compositions P = 0.1 (dotted line), 0.5 (short-dashed line), 0.9 
(long-dashed line). To easily compare the annihilation rate for the Earth to the one for 
the Sun, the plot shown here for the Sun refers to the effective annihilation rate, defined 
as F^ times the ratio {R^/RqY = 1-8 x 10~^, where is the radius of the Earth and 
Rq is the Sun-Earth distance. 

Figure 8 Differential neutrino spectra dN^/dEi, for neutralino annihilation in the 
Earth, as a function of the neutrino energy Ei,. These spectra are calculated for a neutralino 

composition P = 0.5 and for three different values of neutralino mass: = 20 GeV 
(dashed line), = 40 GeV (solid line) and m,^ = 120 GeV (dash-dotted line). 

Figure 9 Flux $^ of the up-going muons as functions of for annihilation in 
the Earth. The threshold for the muon energy is Ej^ = 2 GeV. In figure (a) the three 
representative neutralino compositions are P = 0.1 (dotted line), 0.5 (short-dashed line), 
0.9 (long-dashed line). In figure (b) P = 0.01 (dotted line), 0.5 (short-dashed line), 0.99 
(long-dashed line). The horizontal line in figure (a) represents the Kamiokande upper 
bound 4.0 • lO-^^cm-^s-i (90 % C.L.) [2]. 

Figure 10 Differential muon spectra dNn/dE^ as a function of the muon energy E^j,. 
The solid histogram is the muon spectrum obtained with our MC simulation for P = 0.5 
and = 120 GeV. The solid line represents the muon spectrum calculated in the no- 
straggling approximation. The dash-dotted histogram is the muon energy distribution for 
the background of atmospheric neutrinos. The spectra refer to muons entering the detector 
within a cone of half-aperture of 30° centered at the nadir. 



12 



Figure 11 Muon angular distribution dN^/d cos 9 as a function of cos^, where 9 is 
the nadir angle. The solid histogram is the angular distribution of the signal obtained with 
our MC simulation for P = 0.5 and = 120 GeV. The dashed histogram is the muon 
angular distribution for the background of atmospheric neutrinos. 

Figure 12 Exposure necessary to have a 4a effect with a signal of at least 4 events, as 
a function of neutralino mass . The figure refers to a signal coming from the Earth. The 
three representative neutralino compositions are P = 0.1 (dotted line), 0.5 (short-dashed 
line), 0.9 (long-dashed line). In figure (a) Atmin is displayed keeping the half-aperture 
angle of the detector fixed at 30° from the center of the Earth. Figure (b) displays the 
exposure Atopt obtained by employing the optimization procedure. 

Figure 13 Flux $^ of the up-going muons as functions of for x~X annihilation in 
the Sun. The flux is orthogonal to the Sun direction and the threshold for the muon energy 
is E^^ = 2 GeV. The three representative neutralino compositions are P = 0.1 (dotted 
line), 0.5 (short-dashed line), 0.9 (long-dashed line). The horizontal line represents the 
Kamiokande upper bound 6.6 • lO-^^cm-^s-^ (90 % C.L.) [2]. 

Figure 14 Effective exposure {eAt) necessary to have a 4a" effect with a signal of at 
least of 4 events, as a function of neutralino mass m^. The figure refers to a signal coming 
from the Sun. The three representative neutralino compositions are P = 0.1 (dotted line), 
0.5 (short-dashed line), 0.9 (long-dashed line). The half-aperture of the detector is fixed 
at 5° from the Sun direction. 
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